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Stability Analysis of Optimal Velocity Model for Traffic and Granular 
Flow under Open Boundary Condition 

Namiko MitarajQ and Hiizu NakanishiQ 

Department of Physics, Kyushu University 33, Fukuoka 812-8581 

We analyzed the stability of the uniform flow solution in the optimal velocity model for traffic 
and granular flow under the open boundary condition, ft was demonstrated that, even within the 
linearly unstable region, there is a parameter region where the uniform solution is stable against 
a localized perturbation. We also found an oscillatory solution in the linearly unstable region 
and its period is not commensurate with the periodicity of the car index space. The oscillatory 
solution has some features in common with the synchronized flow observed in real traffic. 
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The traffic_fk)isz^on a highwayS'i'BBi-' and the granular 
pipe ffowaEHa'ErEH' are similar in a number of character- 
istics. Both consist of discrete elements following dis- 
sipative dynamics, both are quasi one-dimensional sys- 
tems, and spontaneous density waves with the power 
law occur in both systems. To understand the common 
structure of such phenomena, both discrete and contin- 
uous models have been suggested and have succeeded 
in renrn d 1 1 ci n p- worn n a spects of the spatiotemporal pat- 
terns .iBBBBSBB 

In recent years, a car-following model, termed the-L'on- 
timal velocity model" (OV model), was proposedB't3 
The model is realistic enough to be able to reproduce 
a spontaneous traffic jam as in real traffic, but simple 
enough to be regarded as a model of granular pipe flow 
without much modification. In fact, a model similar to 
the OV model has been proposed as a granular flow 
model and it has been demonstrated to exhibit power 
law behavior in the power spectra of the density wave in 
certain situations .eI* 

Even though a number of analyses of the OV model 



pted un- 
M The 



have been performed, most of them were cc 
der the periodic boundary condition (PBC).Ei 
PBC corresponds to the situation that cars run in a cir- 
cuit, but it is not suitable for most situations in real 
traffic observations. When we consider the situation of 
granular pipe flow, the PBC is quite unrealistic. Fur- 
thermore, from the experimental results of the granular 
pipe flow, it has been suggested that_the condition of the 
downstream boundary is important .0' 

In the present work, we investigate the OV model un- 
der the open boundary condition (OBC), and examine 
the effects of the boundary condition on the stability of 
the system. 

The OV model is a dynamical model for a one- 
dimensional system which consists of discrete elements; 
the elements can be cars or grain particles, but they tend 
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Fig. 1. The optimal velocity model. 



to run at the speed determined by a local configuration. 
When the {n + l)th car precedes the nth car (Fig. 1) 
and the position of the nth car at the time t is denoted 
by Xn {t) , then the dynamics is governed by the equation 
of motion 



Xn{t)^a[U{bn{t))-Xn{t)], 



(1) 



where &„(t) represents the headway of the nth car defined 
by 



bn{t) = Xn+l{t) - Xn{t), 



(2) 



and the dot denotes the time derivative. Here, a is a 
sensitivity constant and U{h) is the optimal velocity that 
the drivers prefer when the headway is b. For C/(6), we 
employ 



U{b) = tanh(fe - 2) + tanh(2), 

as in most other works on the OV model. 
Equation (|l|) has a uniform solution 

Xn{t) —bn + U {b)t + const.. 



(3) 



(4) 



which represents that all of the cars run with the same 
headwajLij^d the same optimal velocity U(b). It is easy 
to showtlUlS that the solution is linearly unstable when 



a < 2U'{b), 



(5) 



where the prime represents a derivative. The behavior of 
the model in the lineaxly unstable region has been stud- 
ied and demonstratedtilBB that the system eventually 
segregates into two regions, namely, the jammed region 
and the free-flow region. The headways in the jammed 
and the free-flow region are determined by a for a given 
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OV function 

For the purpose of analyzing the system behavior in 
the situation where the upper and lower streams are dis- 
tinguished, we employ the OBC, and define it as follows: 
at cc = 0, the cars enter the system at a constant time 
interval b/U{b) and velocity U(b), and the car that is 
farthest ahead, which has the largest index and no car 
to follow within the system, travels as if its headway is 
6, namely, it follows the equation of motion as 

Xf^a[U{b)-Xf], (6) 

until it leaves the system at x = L. 

In order to examine the stability, we perturb the sys- 
tem at t = locally in space and time by shifting the 
velocity of the car with n = at x = L/2; in the ac- 
tual simulations, the system is prepared with the initial 
condition as 

XniO) = bn + L/2, Xn{0) = U(b) for n = ±1, ±2, 

(7) 

xq{0) = L/2, io(0) = C/(6) + e, (8) 
with a small value of e. 
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Fig. 2. Spatiotemporal diagram of the OV model with a = 1.4, 
b = 2, L = 200 and e = 0.1. The system is perturbed at the point 
indicated by the arrow on the bottom line. Positions of cars 
are indicated by dots, whose size represents the headway of the 
car; the large, medium-sized, and small dots represent cars with 
headways smaller than, equal to, and larger than 5, respectively. 
As a result, the original uniform flow region appears shaded, and 
the region with bn > b {bn < b) appears lighter (darker). The 
disturbance recedes towards the upper stream. 



Figures 2 and 3 show two typical spatiotemporal dia- 
grams, where the original uniform region with the con- 
stant headway bn — b is denoted by the shaded region, 
and the region with bn > b and that with bn < b are 
shown by the lighter and the darker region, respectively. 
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Fig. 3. Spatiotemporal diagram with a = 1, 6 = 2, L = 200 and 
e = 0.1. The disturbance spreads in both directions. 



The perturbation is given ai x = L/2 and t — 0. For 
both cases, the parameters are in the linearly unstable 
region, thus, the largest response to the perturbation 
over the entire system grows as time passes. However, in 
the case of Fig. 2 with the larger value of a, the effects 
of the perturbation travel only towards the upper stream 
and eventually exit the system, then, the entire system 
recovers the uniform configuration. On the other hand, 
in the case of Fig. 3 with the smaller value of a, the 
effects of the perturbation extend in both directions and 
the uniform region is eliminated completely. Therefore, 
within the linearly unstable region in the b — a plane, we 
can draw a line that separates the two distinct parame- 
ter regions: the locally stable region where a small local 
disturbance in the uniform state recedes, and the locally 
unstfible region where the local disturbance extends (Fig. 

4).y' 

What is peculiar to these spatiotemporal diagrams is 
that there appears a region with a regular striping pat- 
tern of lighter and darker regions behind the uniform 
flow region. The behavior of the system in this region 
can be seen more clearly in the snapshots of the head- 
way of the cars. In Fig. 5(a), the right end region at 
n > —300 represents the uniform flow, and the fluctu- 
ating region around the part near n w —450 represents 
an alternating sequence of jammed and free-flow regions; 
the upper and lower limiting values of 5„ in this region 
are close to those of the jammed and free-flow regions 
in the PBC. Between the uniform and the alternating 
region, there exists the oscillatory flow region which has 
a smaller amplitude of fluctuation than that of the al- 
ternating region (Fig. 5(b)); this oscillatory flow region 
corresponds to the regular stripes in Fig. 3. 

This oscillatory flow has some interesting characteris- 
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Fig. 4. Phase diagram in the 5 — a plane. The solid line represents 
the linear stability limit above which the uniform solution is 
linearly stable. The solid boxes represent the local stability limit 
determined by Vo = 0, and the open circles denote the boundary 
where Vi = 0. 



tics, as can be seen from Fig. 5: (i) The oscillatory flow 
region appears directly behind the uniform flow, and the 
oscillatory flow is followed by the alternating region of 
the free-flow and the jammed region when the system 
size L is sufficiently large, (ii) The amplitude of the os- 
cillatory flow (the difference between the minimum and 
maximum value of 5„ in the region) is smaller than the 
difference in the headway between the jammed and the 
free-flow regions, (iii) The oscillatory flow region typi- 
cally extends over a hundred cars and persists, (iv) The 
velocity and the headway of the cars oscillate periodically 
around that of the uniform flow, but their periodicity is 
incommensurate with the period of the car index space. 

Because the oscillatory flow appears spontaneously be- 
ing triggered by a local disturbance, it is expected that 
there exists, at least, a metastable solution where the os- 
cillatory flow extends over the entire system. In order to 
find this oscillatory solution, we assume it is represented 
as 



b + J{n + ct), 



(9) 



with the phase speed c; b is the mean headway, therefore 
the function f{n + ct) should have zero mean. Substitut- 
ing eq. (||) into eq. (|l|), we obtain 

c^f'{z)^a[Uib + f{z + l))-U{b+f{z))-cr{z)]. (10) 

This can be solved numerically and is found to have an 
oscillatory solution for a finite range of the phase speed 
c; e.g., c € (0.56,0.64) for a = 1 and 6 = 2. If we 
set the phase speed at the value obtained by the direct 
simulation of eq. (|l|), namely, c = = 0.61 for a = 1 
and 6 = 2, the resulting f{z) coincides with the results 
of the simulation, as is shown in Fig. 5(c). This clearly 
confirms our conjecture that there exists an oscillatory 
fiow solution. 

The configuration seen in Fig. 5(c) can be regarded as 
the solution that connects the two domains, namely, the 
right one filled with the uniform flow solution, and the 
left one with the oscillatory flow solution. The transition 
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Fig. 5. Snapshots of car configuration represented by the b„ vs 
n plot at t = 988 with a = l,b = 2, L = 10000 (the open circles 
connected by the dashed line), (a) The oscillatory flow region 
appears around —360 < n < —300 followed by the alternating 
jammed and free-flow regions around —550 < n < —400. (b) 
The magnification of (a). The oscillatory behavior can be seen 
clearly, (c) The same data of 6„ are plotted with the solid line 
which represents the result from eq. (fUj) with c = Cs = 0.61. 



region, or the domain wall, is moving at speed Vq. Then, 
the boundary in the 6— a plane between the locally stable 
and the locally unstable region is understood as the line 
where the speed Vq of the domain wall connecting the 
two solutions becomes zero; Vb < {Vq > 0) for the 
locally stable (unstable) parameter region. 

Near the point in the 6 — a plane where the local stabil- 
ity limit, defined by Vq = Oi reaches the linear stability 
limit, a new situation arises (Fig. 6). The free-fiow re- 
gion with the lower car density, which emerges at the 
downstream end of the alternating region, invades the 
oscillatory, and eventually, the uniform flow. The veloc- 
ity Vi by which the boundary between the oscillatory flow 
and the low-density free-flow advances is determined by 
the continuity equation of the flux as 



bi 



U{bi) - Vi U{b) - Vi • 



(11) 



where 6; is the headway of the low-density flow region. 
Therefore, when V/ > 0, the uniform solution is invaded 
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Fig. 6. Spatiotemporal diagram with a = 1.3, b = 2.5, L = 1000 
and e = 0.1. The large headway region is seen to be invading 
the oscillatory flow region. 



by the low-density flow, even when Vq < 0. Vi depends 
on bi, and 6; does not appear to be determined solely 
by a. It has, however, been numerically shown that 5/ 
is always smaller than bf, which is the headway of the 
free-flow under the PBC. Thus, eq. ( [Tl| ) with bi ~ bf 
yields the stability limit, which is also shown in Fig. 4. 

In summary, we have demonstrated that the uniform 
solution of the OV model can be stable against a local 
perturbation, even in the linearly unstable region if we 
employ the OBC. We have also demonstrated that there 
exists an oscillatory flow solution in the linearly unstable 
region and the local stability limit is determined by the 
speed of the domain wall connecting the uniform flow 
solution and the oscillatory flow solution. 

Before concluding, let us make a few comments regard- 
ing the oscillatory solution. This is a new solution, which 
has not been noticed in the weak nonlinear analysis near 
the linear stability limit, where eq. (|^) is reduced to the 
Korteweg-de Vries (KdV) or modified KdV equation.E3 
The system chooses a particular phase speed Cs in spite 
of the finite continuous band of the phase speed being 
allowed. This oscillatory solution has some characteris- 
tics in common with synchronize! 
observed recently in real trafficp'l 
tent fluctuating flow excited by a localized perturbation 
in a high density region where the jammed phase can be 
nucleated. 
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